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The Josephson effect describes the generic appearance of a supercurrent 
in a weak link between two superconductors. Its exact physical nature how¬ 
ever deeply influences the properties of the supercurrent. Detailed studies of 
Josephson junctions can reveal microscopic properties of the superconduct¬ 
ing pairing (spin-triplet correlations, d-wave symmetry) or of the electronic 
transport (quantum dot, ballistic channels). In recent years, considerable 
efforts have focused on the coupling of superconductors to topological insula¬ 
tors, in which transport is mediated by topologically protected Dirac surface 
states with helical spin polarization (while the bulk remains insulating). Here, 
the proximity of a superconductor is predicted to give rise to unconventional 
induced p-wave superconductivity, with a doublet of topologically protected 
gapless Andreev bound states, whose energies varies 47r-periodically with the 
superconducting phase difference across the junction. In this article, we re¬ 
port the observation of an anomalous response to rf irradiation in a Josephson 
junction with a weak link of the 3D topological insulator HgTe. The response 
is understood as due to a Jyr-periodic contribution to the supercurrent, and its 
amplitude is compatible with the expected contribution of a gapless Andreev 
doublet. 

The helical nature of the topological surface states, where the spin is locked perpendic¬ 
ular to the momentum^, is predicted to give rise to exotic superconductivity when coupled 
to the conventional pairing potential of a s-type superconductor. The broken spin rota¬ 
tion symmetry allows the appearance of triplet p-wave correlations and of gapless Andreev 
bound states, regardless of the microscopic details of the theoretical modeP®. In a 3D 
topological insulator (TI) based Josephson junction (JJ), in which superconductivity is 
induced by the proximity effect of a s-wave superconductor (see Figj^), Andreev bound 
states appear in the induced gap A* that can be pictured as (see Fig|^) a single topo¬ 
logical Andreev doublet (depicted in blue) that occurs at transverse momentum ky = 0 
and is immune to back-scattering (thus has perfect transmission), and non-topological 


2 


oblique modes {ky ^ 0, depicted in red) that are expected to have lower transmissions. 
The topological protection of the zero mode constitutes a superconducting analogue to 
Klein tunneling. As depicted in the peculiarity of this topological doublet is its 

dvr-periodicity (or equivalently a contribution sin 0/2 to the supercurrent) with respect 
to the superconducting phase difference 0 across the junctioi#*^. The non-ambiguous 
observation of such gapless states is regarded as an important experimental signature of 
the unconventional superconductivity in topological insulators, but no robust evidence 
has been reported yelP®^. A major hindrance could be the coexistence of residual bulk 
conductance or of a large number of gapped conventional modeJ^^Ull Furthermore the 
finite lifetime of the positive energy branch prevents observation of a dvr-periodic Joseph- 
son effect in stationary measurements as various relaxation mechanisms can restore a 
2vr-periodicity for the current-phase relation (CPR)^^^. The dvr-periodic Josephson ef¬ 
fect can therefore be unveiled more easily by the dynamics of the junction. To reveal the 
periodicity of the Josephson supercurrent, an rf driving current A/ is added to the dc drive 
to induce the so-called Shapiro stepJ^Sl. When the dynamics of a conventional JJ is phase- 
locked to the rf drive, steps of constant voltage appear in the l-V characteristic of the 
junction for voltages 14 = nhf j2e where n G Z is the step index. However, for a purely 
dvr-periodic supercurrent, only a sequence of even steps should be observed. In the case of 
nanowire^, signs of the disappearance of the first step {n — 1) have been reported and 
attributed to a theoretically expected topological phase transition driven by a magnetic 
field along the axis of the nanowire, although the topological state in the nanowires has 
yet to be identified in the normal transport regime. Even though it is known that other 
systems (JJ between p/d-wave superconductorJ^^*^ could lead to similar effects without 
the need for topological protection, such effects have yet to be observed. In this article, 
we study HgTe, a genuine 3D topological insulator whose topological properties have been 
established independentljPIES, and observe an anomalous doubled Shapiro step appearing 
at low frequency (equivalently a missing n = 1 step). While several other mechanisms 
(non-linearities, capacitance effects, higher harmonics in the current phase relatiorP^Edl) 
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are known to cause the appearance of additional subharmonic steps in the Shapiro re¬ 
sponse, to our knowledge, only the existence of a 47r-periodic contribution / 4 ^ sin 0/2 in 
the total supercurrent can be responsible for the disappearance of odd stepd^. 

Our devices are fabricated from coherently strained undoped HgTe layers of 65 to 
90 nm thickness, epitaxially grown on a CdTe substrate. The band inversion of HgTe 
enforces the existence of topological surface states, while strain opens a gap (~ 22 meV) 
in the bulk of the materiaP^. Previous work has highlighted the high quality of the 
topological states in this materiaP^IES. Quantized Hall plateaus are routinely observed, 
which demonstrate that transport occurs exclusively through the surface states, without 
any detectable parallel conductance from the bulk. The mobility and charge density, 
relevant for our experiments, are evaluated from a Hall-bar produced separately from 
the same wafer as the junctions, and yield typically fi — 1 — 3 x 10^ cm^ s“^, and 
Ue = 3 — 7 X 10^^ cm“^. From these values, we extract a mean free path of I ~ 200 nm. 
The JJs are fabricated by depositing niobium contacts at the surface of a HgTe mesa, 
using standard sputtering and lift-off techniques (see Supplementary information). The 
geometry is shown in Fig{T]A.. Each superconducting contact has a width of 1 to 4 pm, the 
HgTe weak link has a width of kF = 2 pm (corresponding to the width of the mesa stripe) 
and a variable length L ranging from 150 nm to 600 nm. From the electron density Ue, we 
evaluate the number of transport modes N — ~ 100. The niobium of the contacts 

has a critical temperature of Tc ~ 8 K, slightly lower than that of bulk Nb (9.2 K). 
A typical I-V curve obtained at 30 mK is presented in Figj^) and exhibits hysteresis, 
as commonly reportecP^. We find that the critical current of devices with the same 
dimensions varies by about 30%, which underlines the reproducibility and quality of the 
fabrication process. A recurring feature in all devices is the presence of an excess current 
in the I-V curve (see Figj^). For high biases, the I-V curves become linear with an 
asymptote which does not go through the origin but is shifted towards higher currents. 
This excess current is understood as due to the fact that electrons in an energy window 
near the superconducting gap carry twice as much current due to Andreev refiection d^^*^^ l 
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It thus illustrates the presence of Andreev reflections at both S-TI interfaces. Such an 
excess current has been previously observed for superconducting point-contact#^*^, but 
is not commonly reported in thin film structures presumably due to the presence of 
elastic scattering. This emphasizes the high quality and reproducibility of our devices in 
agreement with our previous observation#*^®^. 

We now turn to the study of the response of these devices to rf irradiation and high¬ 
light the existence of a dvr-periodic supercurrent. To this end, we focus on three devices 
produced from the same wafer, for which the width of the junction is set to W = 2 pm, 
for nominal lengths of L = 150nm, 400nm and 600 nm (see Supplementary information). 
The experiment described below has been repeated on more than ten devices, made out of 
three different wafers with similar characteristics, in three different measurement setups, 
all yielded similar results (see Supplementary information). 

In order to observe the Shapiro steps, the sample is irradiated with a radio-frequency 
excitation via a coaxial line, the open end of which is adjusted to be around 1 mm 
from the sample. In this geometry, frequencies in the range of 2 to 12 GHz are easily 
accessible, but the rf power supplied to the sample is not calibrated. Under rf irradiation, 
we observe the appearance of Shapiro steps in the I-V characteristic at quantized voltages 
Vn = nhf /2e, where n e Z is the step indexP*. In contrast to the standard Josephson 
junction response, with steps at each n, we And at lower frequency that the n = 1 step 
is missing. To illustrate this anomalous Shapiro response of our junctions, we present 
three I-V curves corresponding to three different excitation frequencies in Figj^ (for the 
junction with L = 150nm). The applied rf power is chosen such that all curves display 
similar critical currents, the full range of rf power will be discussed later. For a high 
frequency / = 11.2 GHz, one typical I-V curve is plotted as a blue line (with voltage 
normalized to hf /2e). Several steps are clearly visible with step height hf /2e. At lower 
frequencies / = 5.3GHz (green line), higher order steps are visible but a clear reduction 
of the amplitude of the n = 1 step occurs. For a frequency of / = 2.7 GHz (red line), this 
first odd step is fully suppressed, showing an anomalous first step at hf /e. The presence 
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or absence of the n = 1 can be conveniently detected by binning the measurement data 
according to the voltage (with a 0.25|^ bin size). The resulting histograms of the voltage 
V are presented as bar plots in Figj^. For Vn = nhf /2e with n integer, Shapiro steps 
appear as peaks in the bin counts, the amplitude of which then reflects the length of 
the current step (in nA). For / = 11.2GHz (right graph), all steps emerge clearly from 
the background. On the left side, for / = 2.7 GHz, the peak at V = hf/2e is absent, 
reflecting the suppression of the n = 1 Shapiro step. This anomalous behaviour of the 
Shapiro steps constitutes the main finding of this article. Below, we carefully analyze its 
origin and conclude that it indicates the existence of a 47r-periodic contribution to the 
supercurrent. 

We now examine the crossover from high to low frequency, for which the first odd 
Shapiro step n — 1 progressively disappears. To this end, we scan the presence of Shapiro 
steps for a range of rf powers at fixed frequencies and generate two-dimensional color plots 
of the bin counts at the voltage V (which indicates the current-height of the Shapiro step 
when present) as a function of the voltage V and rf current A/. As shown in Figj^ (for 
the junction with L = 150 nm), such plots reveal the presence of Shapiro steps as maxima 
at constant quantized voltages (horizontal lines). Let us first examine measurements 
taken at / = 11.2GHz. (Figj^). At A/ = 0, a single maximum at G = 0 reflects the 
presence of a supercurrent. As A/ increases, Shapiro steps progressively appear, starting 
from low values of n, while the amplitude of the supercurrent (n = 0) decreases and 
eventually vanishes. At higher powers, the steps show an oscillatory pattern, reminiscent 
of Bessel functions occurring in the voltage bias cas (j^^ * ^^ l Horizontal linecuts at constant 
voltages give access to the amplitude of the first steps {n = 0,1, 2,3,4), presented in 
the lower panels of Figj^ as a function of rf current A/- For high frequencies such as 
/ = 11.2 GHz, our device exhibits the conventional behavior which is seen in various other 
systems (carbon nanotubed^, graphen^, or Bi 2 Se^weak links), that always (regardless 
of frequency) show a clear presence of the n = 1 step. The case of atomic contacts (with a 
few ballistic highly-transparent modes) is particularly well understood, and also exhibits 


6 




a strong n = 1 Shapiro resonance in excellent agreement with theoretical modeld^^®. In 
the Supplementary Information, we provide additional measurements on graphene-based 
devices (another example of 2D Dirac material), that also shows this standard behavior 
in all accessible regimes. 

In contrast to the conventional Shapiro features commonly reported, our HgTe-based 
junctions exhibit a very clear vanishing of the hrst step n = 1 when the excitation fre¬ 
quency / is decreased. Measurements at / = 5.3 GHz show that the first step is suppressed 
below a certain value of 7^/ (indicated by the red arrow), and that it is completely ab¬ 
sent at / = 2.7 GHz. In the oscillatory regime at higher rf currents, a suppressed first 
oscillation (dark fringe indicated by the dark grey arrow) becomes clearly visible at low 
frequency, demonstrating the range of influence of the vanishing first step on the rest of 
the pattern. In the lower panels, a complete suppression of the first step or disturbances 
in the oscillations at higher rf currents can similarly be observed. This crossover has been 
observed on all working devices, up to 800 mK, which is the highest stable temperature 
accessible in our fridge. In some cases, hysteretic behavior at low temperatures hinders 
the observation of low-index steps (see Supplementary Information). However, impor¬ 
tantly, biasing instabilities and sudden current switches (such as the ones observed in the 
hysteretic regime) can be excluded as a mechanism for the missing n = 1 step, as the 
same features are seen in measurements of a junction in which bistability is suppressed 
by a shunt resistor (see Supplementary Information). 

In opposition to a missing n — 1 step, additional subharmonic steps (for n = p/q 
fractional value) are often observecP^® as a consequence of non-linearities, capacitance 
effects or higher harmonics in the GPR. Such higher harmonics have been predictecP and 
detected^ in our junctions. At higher frequencies, we indeed observe half-integer steps 
{n = l/2,3/2,..., see Supplementary Information) but they clearly appear in a different 
regime from where we observe the missing n = 1 step. 

The presence of a dvr-contribution in the supercurrent I 4 ^,-sine/)/2 is the only known 
mechanism to result in the observed doubling of the Shapiro step size. As already men- 
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tioned, microscopic models based on Bogoliubov-de Gennes equations have predicted such 
a dTT-periodic contribution in the CPRp^, which originates from the presence of a gap¬ 
less topological Andreev doublet. This anomalous CPR can then be supplemented with 
the Josephson equation on the time-evolution of the phase difference to simulate the 
dynamics of such a system. This dynamics is captured in the extended RSJ model of 
Dominguez et al^. It takes into account the presence of a sin 0/2 contribution in the 
supercurrent and explains the crossover between the two frequency regimes by the highly 
non-linear dynamics of the junction. When a small dyr-periodic contribution 14 ^ sin 0/2 
is superposed on a large 27r-periodic supercurrent I 2 tt sin 0 in the current-phase relation, 
the latter dominates the high-frequency Shapiro response, but the weak dvr contribution 
is revealed at low frequencies by doubled Shapiro steps (see Supplementary Information). 
Doubled Shapiro steps are observed only when the driving frequency / becomes smaller 
than the characteristic frequency 04 ^ = (with the normal state resistance of the 

device). This frequency scale based on the amplitude of the dvr supercurrent is expected 
to be much smaller than the typical Josephson frequency scale fj = ~ 53 qHz 

for the 150 nm long junction), as <C hn — Ic- In order to estimate we introduce 
two indicators Q 12 and Qu as follows. From the maximum amplitude of the first lobe of 
each step, denoted by Wn, n e Z, (see Figj^ where the measurement is indicated for the 
n = A step), we define and compute the ratios Q 12 = Wi/w 2 i <534 = vjs/w 4 , and plot them 
as a function of the rf excitation frequency (Figj^. Despite some scattering, we observe 
a clear decrease of Q 12 towards 0 with decreasing frequency, while Q 34 remains constant 
around 1, for all lengths. For the shortest junction (150 nm) Q 12 reaches a value of 0.05 
around 2 GHz, and the first step n = 1 is invisible. For comparison, we have also plotted 
the boundaries (grey dashed lines) between which the ratios Q 12 and Q 34 vary in the stan¬ 
dard RSJ modeP^*^ (with only a sin 0 component in the supercurrent, see Supplementary 
Information). While the ratio Qs 4 remains close to grey region, the behavior of Q 12 is not 
properly described. Assuming the validity of the above criterion, one can evaluate the 
number of dvr-periodic channels. We estimate f 4 -,^ = 4.5 — 5 GHz and = 250 — 300 nA 






for the 150 nm junction, and / 47 r = 4 GHz and = 50 — 70 nA for the longer junctions 
(400 and 600 nm). One can compare these valnes with the maximum supercurrent carried 
by one channeP^, given by eAj/^ per channel where Aj can be estimated from the decay 
of Ic with temperature (see Supplementary information). With A, = 0.35 meV (150 nm) 
and Aj = 0.1 —0.15 meV (400 and 600 nm), we estimate that the dvr-contribution amounts 
to that of 1-3 channels which is compatible with the presence of one topological mode in 
our system, despite uncertainties on the exact valne of / 47 r and A*. Results on three other 
jnnctions with different parameters have been compiled (see Supplementary Information) 
and are consistent with this estimate. 

Finally, one might also suspect that the dvr-periodic contribution stems from Landau- 
Zener transitions occnrring at the anticrossing (for 0 = vr[2vr]), causing some highly- 
transparent 2vr-periodic states to behave effectively as 4vr-periodic, in the absence of truly 
4vr-periodic modes. In a single mode modeP^, one can numerically show that the quanti¬ 
zation of the Shapiro steps is lost when the Landau-Zener tunneling probability is lower 
than 1: the Shapiro steps split in two and depart from their quantized valnes 14, even), 
and eventnally disappear for probabilities below 0.7. We do not experimentally observe 
snch effects in any accessible regime. Assuming the validity of this specific model, an 
upper bound on the possible energy splitting 25 between positive and negative energy 
branches can be evaluated from the Landau-Zener transition probability. Given our ex¬ 
perimental resolution, we obtain the upper bound 5 < 9peV for the 400 and 600 nm 
jnnctions, and 5 < ISpeV for the 150 nm one (see Supplementary Information). This is 
mnch smaller than the energy scale given by the temperature (70peV at 800 mK) and 
corresponds in both cases to a transmission > 0.994. Besides, there are no reports of 
missing odd Shapiro steps due to Landau-Zener transitions in highly ballistic junctions 
reported to date. 

Interestingly, only the first step n = 1 is missin^^, and not the following odd steps 
n = 3,5,.... This feature is not well understood, but one possible explanation is given 
by enhanced relaxation (due for example to coupling to the continuum of states above 
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the superconducting gap), with a characteristic time scale that decreases as voltage 
increase d^^ * ^^ ! However, most models assume a voltage bias of the junction and a more 
detailed analysis of the current bias case is needed. 

To conclude, we have presented robust evidence for a dvr-periodic contribution to the 
supercurrent flowing in Josephson junctions based on the 3D topological insulator HgTe. 
The consistency of the measurements in Hall bars and JJs signals that our devices are 
well-controlled, with well-defined proximity-induced superconducting HgTe contacts con¬ 
nected via a ballistic HgTe surface. Under rf irradiation, a suppression of the first Shapiro 
step is observed at low frequencies and low magnetic fields, for a wide range of temper¬ 
atures (up to 800 mK), which we attribute to the existence of a dvr-periodic component 
in the supercurrent. The study of its order of magnitude and of Landau-Zener transi¬ 
tions reveal that these experimental observations are compatible with the presence of a 
few dvr-periodic gapless Andreev bound states. Such states would likely stem from the 
topologically protected gapless Andreev bound states, but could also originate from triv¬ 
ial ballistic states. Further investigations are required to conclusively demonstrate the 
relationship of these observations to Majorana physic^®. Besides, these observations 
with a 3D TI of strained HgTe are very encouraging for future experiments in which the 
weak link would consist of narrow HgTe quantum wells that exhibit the quantum spin 
Hall effecPl, in which the total number of transport modes should be reduced to a few. 
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FIG. 1. Geometry of the Josephson junction and predicted Andreev spectrum - 

(A) Artist view of the Josephson junction. Mesa stripes of HgTe (represented in mauve) are 
patterned on the CdTe substrate (blue), with a width W = 2pm. Nb contacts (in yellow) 
are added at the surface, with a width of of 1 to 4 pm, separated by a variable distance L. 

(B) Typical energy spectra £((/>) of the Andreev bound states in a 3D TI based junction, as a 
function of the phase difference (j) in the JJ. In blue is depicted the gapless dvr-periodic topological 
mode, corresponding to transverse momentum ky = 0 and contributing to the dvr-supercurrent 
/ 47 r sin(/)/2. In red, gapped modes correspond to ky ^ 0 and contribute to the 2vr-supercurrent 
I 27 T (containing sin^ and higher order harmonics). 
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FIG. 2. I-V curves of the L = 150 nm junction) - (A) I-V curve in the absence of rf 
irradiation for the two sweep directions, taken at base temperature T ^ 30 mK. The asymptotes 
(grey solid lines) do not cross the origin, emphasizing the presence of an excess current. (Inset) 
Detailed view of the I-V curve, that exhibits hysteresis between the upward and downward 
sweep direction. (B) Shapiro steps for three different frequencies measured at T 800 mK. The 
plotted voltage scale is in normalized units hf /2e to highlight the formation of Shapiro steps 
in the I-V curve in the presence of rf irradiation. For a high frequency / = 11.2 GHz (blue 
line), all steps are clearly visible for voltages Vn = '^^5 n ^ Z (up to |n| > 12, but only the 
first three are shown for the sake of clarity). For an intermediate frequency (/ = 5.3 GHz, blue 
line), the first step (n = 1) is noticeably reduced. At low frequency (/ = 2.7GHz, red line), 
the first step is fully suppressed, while all other steps remain visible. (C) Bar plots obtained 
by binning the measurement data according to voltage, for / = 2.7 GHz and 11.2 GHz. The 
Shapiro steps appear as peaks in the bin counts for Vn = n G Z. While all steps are visible 
for / = 11.2 GHz (right graph), the first Shapiro step {n — 1) is absent at / = 2.7 GHz (left 
graph). 
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FIG. 3. 2D plots of the bin counts and Shapiro step amplitudes for the L = 150 nm 
junction) - Differential conductance is plotted as a function of the voltage V (normalized 
by the step size at frequency / hf /2e) and rf current I^f for rf excitation frequencies / = 
2.7GHz,5.3GHz, 11.2GHz in (A), (B) and (C) respectively, (upper panels) Shapiro steps are 
identified as maxima for constant voltages Vn (white dashed lines emphasize n = 0,1,2). For 
/ = 11.2 GHz, all steps are visible. When frequency is lowered (/ = 5.3 GHz), the first odd step 
{n = 1) is absent up to a rf excitation indicated by the red arrow. Finally, at / = 2.7 GHz, 
the first step is completely invisible up to the crossing point which marks the beginning of the 
oscillatory regime at high rf currents. A dark fringe (indicated by a dark grey arrow) is observed 
at finite voltages in the oscillating pattern concomitant with the missing n = 1 step, (lower 
panels) Horizontal line-cuts through the previous colormaps give access to the amplitudes of 
steps 0 to 4. While all Shapiro steps are clearly visible at high frequencies, the step n = 1 
progressively disappears as / decreases. From these plots, we access the maximum widths Wn 
of each step (see the example of at f = 2.I4GHz). For clarity, the different curves are offset 
by 0.4, 0.6 and 1.7pA for (A), (B), (C) respectively. 












































C-f- 


CO 


FIG. 4. Ratios of step widths Qu and Qs4 vs frequency / - For each length L of the 
JJ, we calculate the ratios of step amplitudes Qu = wilw 2 and Qs 4 = wslw 4 ^ and plot them 
as a function of the rf frequency. Qu shows a very clear decrease as frequency / is lowered. A 
minimum around 0.05 is obtained for the 150 nm junction, but we observe that this minimum 
tends to increase with the length L of the junction. In contrast, even if the measurements 
show some scattering, the ratio of higher order steps Qs 4 does not show significant variation. 
For comparisons, we evaluated Qu and Q 34 from two conventional RSJ models, and show the 
results as a grey area (see Supplementary Information). 
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47r-periodic Josephson supercurrent in HgTe-based topological 

Josephson junctions 


Supplementary Material 

In this supplementary online text, we first detail the sample preparation 
and give further results on the properties of the Josephson junctions. Then 
additional results are presented, obtained from different HgTe-based devices. 
In particular, the issue of subharmonic steps and hysteresis are discussed. 
In a third part, we present the Shapiro response of graphene-based devices 
measured for comparison with the HgTe samples. Finally, we discuss elements 
of theory on RSJ simulations (with and without Jvr-periodic supercurrent) and 
Landau-Zener transitions. 


I. ADDITIONAL CHARACTERIZATION MEASUREMENTS 

A. Sample preparation and layer characterization 

Bulk HgTe layers are grown by molecular beam epitaxy on a CdTe substrate. Prior 
to the fabrication of the JJs, the transport properties of each layer are characterized by 
the measurement of longitudinal and transverse (Hall) resistance in a Hall bar geometry. 
From the longitudinal resistance at zero magnetic field, one can extract the mobility of 
the layer, while the density is obtained from a linear fit of the Hall resistance between 0 
and 500 mT. The layers being very similar to the ones presented in references |2T] and [22l 
we refer the interested reader to these references where the measurements are discussed 
in detail. 

Using a Ti/Si02 etch mask, the samples are patterned via Ar"*" ion beam milling to obtain 
2 pm wide HgTe stripes. The Nb contacts (60 nm Nb + Al/Ru cap) are sputtered on 
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top of the HgTe layer, after a short Ar+-milling step in order to remove any adsorbant 
or reaction product left on the HgTe mesa after exposure to air. The top-layer of HgTe 
below the contact may be disordered due to the deposition of the Nb contacts. We assume 
that this top-surface is proximitized by interaction with the Nb, but that the HgTe weak 
link between the contacts remains in the ballistic regime as highlighted in the main text. 



FIG. SI. SEM pictures of two typical devices - (A) Colorized SEM pictures, highlighting 
the HgTe mesa (in red) and the Nb electrodes (in blue), in the lateral junction geometry pre¬ 
sented in Fig.l (main text). (B) Close-up on the jnnction itself. Nb sidewalls are visible on the 
electrodes, and are a consequence of sputter deposition on the sidewalls of the resist pattern 
used to dehne contacts. 

In FigJSTI SEM pictures of a typical junction are shown. From these SEM pictures, 
it is observed that the nominal length (used in the article) probably overestimates the 
physical length of the JJs. We evaluate the physical length to be 60 to 70 nm smaller 
than the nominal length. 


B. Temperature dependence 


We present in Figj^ the temperature dependence of four quantities. In FigjS^, the 
critical current A and retrapping current A are shown (in red and blue dots respectively). 
Both decay on a typical scale of 1-2 K. The hysteretic region, characterized by A > A 
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takes place at roughly T < 800 mK in most samples. Beyond this temperature, Ic = Ir 
and no hysteresis is observed in the I — V curve. In FigjS^, the excess current Ig^c and 
the normal state resistance Rn of the jnnction are plotted as a fnnction of temperatnre 
T. First, the resistance Rn (red dots) exhibits a jump (from 31 to 51 fl) indicating the 
superconducting transition of the bulk Nb contacts. These measurements were taken in a 
dilution fridge for which stable temperature control was only possible below 1 K. Higher 
temperature measurements were collected while the system warmed up allowing only a 
few reliable measnrement points snch that the exact transition temperature is not known. 
From other measnrements during cool down, we obtain a critical temperature around 
Tc — 8K, slightly smaller than that of pure high-quality Nb (9.2K) and compatible with 
the measurements presented here. Simultaneously, the excess current is measured and 
exhibits a jump from zero to a finite valne (around 4pA) as the Nb contacts become 
snperconducting. This signals the presence of Andreev reflections at both interfaces of 
the jnnctions for all temperatnres T < Tc. 


The temperature dependence of the junction with L = 150 nm presented in FigJS^can 
be analyzed using the approach developed by Tkachov et alW. Each mode (indexed by the 
transverse component of the wavevector ky) is described by an angle-dependent transmis¬ 
sion To where sind = ky/kp. The distribntion of transmission To reflects Klein tnnneling 
through a barrier (characterized by a parameter Z in a BTK-like approaclP^, with topo¬ 
logical protection of the mode ky = 0, {To=q = 1). This model also takes into acconnt 
induced superconductivity in the HgTe reservoir using a McMillan tunneling approach!^, 
which allows for the evaluation of the induced gap A*. Unfortunately, this perturbative 
approach breaks down for high temperature when the induced gap A* becomes comparable 
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FIG. S2. Temperature dependence - (A) Temperature dependence of the critical and retrap¬ 
ping current Ic and A. The hysteretic region is seen for T < 800 mK typically, for which < Ic- 
(B) Temperature dependence of the excess current lexc and normal state resistance Rn- The 
superconducting phase transition in the Nb contacts is seen here as a shift in the resistance 
and typically takes place at Tc ~ 8K. Simultaneously, the excess current Rxc exhibits a jump 
from zero to a finite value, thus reflecting the presence of Andreev reflections at the interfaces 
of the junction for temperatures T <Tc. 


to the Nb gap A^b (that decreases with temperature). 


rim 
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1 -Tesin"! 
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cos^ 9 ^3 ,, 

-—• 2a/n _L 72V ^ r/Ajvfe, Aj = r(i - 7 + -7 ) 

1 — sin 9/{l + Z^) 2 


( 1 ) 


Eq.Q has three fit parameters: the tunneling strength F (or equivalently the induced 
gap Aj), the barrier parameter Z and the number of modes N. As seen on Fig, S3 A., the 
agreement below 1 K is very good for the L = 150 nm junction, but the fit diverges rapidly 
at high temperature. We obtain the following fitting parameters F = 0.4meV, A, = 
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FIG. S3. Fit of of the critical current vs. temperature - (A) Experimental data for the 
critical current Ic{T) as a function of temperature T is presented as blue + symbols, for the 
150 nm JJ. The solid line is obtained using Eq. 0 for three-dimensional topological Josephson 
junctions, with = 0.35 meV, N — 117, Z = 0.1. (inset) Complete temperature range. Beyond 
IK, data is taken as the fridge warms up, resulting in inaccurate measurements of T. The fit 
departs from experimental data as the perturbative model breaks down. (B) The temperature 
dependence of Ic{T) is presented for the 400 and 600 nm junctions, in red and green respectively. 
Experimental data is showed as + symbols while solid lines present fits with parameters = 
0.13meV, N = 48, Z = 1 (600 nm) and A^ = 0.13meV, N = 33, Z = 1 (400 nm). 


0.35 meV, Z = 0.1, = 117. The number of modes is then in agreement with independent 

estimates from Hall bar measurements, and yields an estimate of A^ that can be used to 


estimate the number of 47r-periodic modes (see section IIB). 


For longer junctions (L = 400 — 600nm), the agreement is not as good as for the 150 
nm long junction presented here. The evaluation of the temperature scale over which 
the decay of Ic is observed gives an estimate of A^ = 0.1 — 0.15meV lower than for the 
junction with L = 150 nm. The amplitude of the critical current Ic decreases rapidly with 
length, and yields a number of modes < 50 which seems unreasonably small given the 
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estimate of the density from a separate Hall-bar, or from the value of N in the 150 nm 
long device on the same sample. It could indicate that these junctions depart from the 
short junction limit I <C ^ where ^ is the coherence length in the system. The natural 
coherence length is typically between 250 and 1000 nm in our systenP^l. For 

systems with mean free path I ^ ^o, the relevant length is in fact ^ = \/^ in the range 
of 250 to 550 nm which is compatible with our findings. 


C. Excess current 

The excess current in our data is indeed strikingly present in all our samples. We 
take it as an additional indication of the systematic reproducibility of our data within a 
large variety of samples of different fabrication runs and measured in different cryostats 
at different locations. The excess current in a Josephson junction is obtained for voltages 
beyond 2A, as discussed by Blonder et alW^, and reaches in principle twice the value for a 
NS contact. However, the amplitude of the excess current depends on elastic scattering at 
the interfaces and as such is not a sufficient measure. The presence of gapless modes can 
in principle be detected by a halved onset on the excess current: the asymptotic regime 
is reached for a bias eV on the order of the gap A and not 2A in the conventional case. 

To analyze more precisely this behavior, we plot in the right panel of FigjS4| the differ¬ 
ence lexciy) = I ~ y/Rni with lexciV) ^ 2Ajv6. We observe two inflexions 

around V ~ 1.4 mV and V ~ 2.4 mV, that could be related to A^b and 2A]sib- However, 
these features are relatively weak, and no clear transition at eV ~ Ajvfe is observed. 

The fact that we cannot observe a clear asymptotic regime at A is however not un¬ 
expected as the previous model is only directly applicable to N-S contacts, in which a 
single superconducting contact is involved and thus the dynamics of Josephson effect is 
not present. In a Josephson S-S contact, the observed currents and voltages below 2A 
are mixed with the time-averaged quantities of the Josephson effect. The I-V curve thus 
reveals a time-averaged part due to the Josephson effect and a static part due to the 
Andreev reflection processes. To get experimental access to the excess current part, it 
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FIG. S4. Excess current and midgap states - A) The I-V curve of the 150 nm junction is 
presented in blue and red line (for the two sweep directions). Two asymptotes are presented as 
a grey line, that do not cross the origin, demonstrating the presence of an excess current. Two 
horizontal dashed line have been added as guidelines and represent possible onsets at eV ~ Ajvfe 
and eV ~ 2Ajvfe. B) Excess current Iexc{V) — I — V/Rn as a function of voltage V. Inflexions 
are visible around V ~ 1.4mV and V ~ 2.4mV. At high bias, Iexc{V) — 5.3pA. 

is necessary to separate these contributions. The contributions due to the Josephson ef¬ 
fect are strongest at V = 0 and at finite voltages close to V = 0 but gradually decay 
in amplitude with increasing voltage. In contrast the excess current is present out to 
high voltages beyond 2A and perhaps beyond A. Careful experimental work is therefore 
needed to disentangle these two processes from the experimental data. Furthermore, the 
theory has been developed for nanowires (single-band and short junction limit) with a 
simple representation of induced superconductivity, by a unique gap A. The fact that our 
system is a 2D surface state rather than a nanowire is expected to have little consequence. 
However, we expect to have a very different density of states due to the presence of two 
superconducting gaps in proximity with each other (the bulk Nb gap Ajvb and the induced 
gap Aj). The literatur d^^B^ I shows for example the presence of features at both Ajvfe and 
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Aj, as well as different densities of states in regions of energy Ai < e < Apfi, and e > Apfi,. 
The onset of the excess cnrrent in onr measnrements is of the same order of magnitude 
as Ajv6 — 1 meV, and the role of the induced gap (evaluated around Aj ~ 0.1 — 0.4 meV) 
remains experimentally not accessible. 


II. ADDITIONAL RESULTS ON HgTe-BASED JUNCTIONS 


A. Frequency dependence and half-integer steps 


In FigjS^ we present datasets measured for a different sample than the one presented 
in the main text with length L = 400 nm. In particular, we show the transition from 
a doubled Shapiro step (at / = 3.34 GHz) to a regime in which the first step is fully 
recovered. At high frequencies (for / > 6.2 GHz), we observe the appearance of new half¬ 
integer steps, at voltages given by the step indexes n = 1/2,3/2,5/2,.... To visulaize more 
clearly the subharmonic steps, we introduce a different way of visualizing our data. The 
differential conductance dl/dV is plotted as a colorscale, as a function of the dc voltage 
V and rf current drive A/. Thus, Shapiro steps appear as maxima [dl/dV diverges) for 
constant voltages, similarly to what is seen in the bin counts presented in the main text. 
In such plots, the information on the step amplitude is lost, but subharmonic steps of 
small amplitudes become more visible. 

At / = 3.34 GHz, the first step (n = 1) is partially suppressed, but is recovered 
as frequency increases. The change in the oscillatory pattern is also clearly visible: at 
low frequency, fast oscillations follow a long region in which step n = 0 (supercurrent) 
is present. As / is increased, oscillations start earlier and with a larger period. These 
features are accounted for by the RSJ model presented in section [TV A[ More interestingly, 
half-integer subharmonic steps (n = 1/2,3/2, 5/2...) become clearly visible on the right 
panels (for / > 6.2GHz). 

As discussed in refl21l it can be reproduced from the RSJ equations by adding a 
capacitive shunt in the circuit (RCSJ modeP^. Another possible mechanism is given 
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FIG. S5. Frequency dependence of the Shapiro response - The differential conductance 
dljdV is plotted as a colorscale, as a function of the dc voltage V and rf current drive I^f for 
different frequencies. The progressive appearance of the initially (partially) missing n — 1 step 
is observed. At high frequencies (/ > 6.2 GHz), subharmonic steps n = become 

24 


visible (see section IIA). 































































by the presence of higher order harmonics in the current-phase relation. As a sin 0/2 
term can suppress the odd steps, a sin 20 term generates subharmonic steps. While the 
geometric capacitance is believed to be negligible (see estimate in reflHI), ^ non-sinusoidal 
current phase relation has clearly been established recently in our junctiond^. 

For the sake of completeness, we present in Figj^a dataset measured at / = 13.2 GHz, 
on the junction with a length L = 150 nm presented in main text. The I-V curve estab¬ 
lishes the presence of two weak half-integer steps at n = 1/2 and 3/2. 



DC current I [/tA] 


FIG. S6. Subharmonic steps - We present an I-V curve obtained under rf irradiation at 
/ = 13.2 GHz. The subharmonic structure (clearly visible in the differential resistance dV/dl in 


Fig,S5) is related to the appearance of subharmonic steps with index n = 1/2, 3/2.... Here in 
particular, one observes the steps n = 1/2 and n — 


B. Summary of measured devices 

In the following table, we summarize the parameters of 6 devices for which a complete 
set of data has been collected. Parameters are as follows : length L, critical current /c. 
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normal state resistance Rn, Josephson frequency fj = , crossover frequency / 47 r, 47r- 

periodic contribution to the supercurrent induced gap A*, maximum current per mode 
io = ^, and finally number of 47r-periodic modes. Starred (*) values indicate devices for 
which the temperature dependence of the critical current has not been measured precisely, 
so that Ai is only estimated from the value of R at base temperature, with large error bars. 
Despite a wide range of parameters explored, all junctions yield an estimated number of 
dvr-periodic modes between 1 and 3. ^ indicates the critical current of the junction with 
L = 150 nm presented in the main text at the time when the Shapiro steps have been 
measured. A previous cool-down of the sample yielded R = 5.6 pA as presented in the 
main text (Fig.2). Aging of the sample has been observed in several samples and could 
explain this discrepancy. 



L (nm) 

Ic (fiA) 

Rn i^) 

fj (GHz) 

hn (GHz) 

hw (nA) 

Ai (meV) 

R (nA) 

N^pi 

A 

150 

3.3^ 

33 

53 

4.5-5 

250-300 

0.35 

90 

3 

B 

400 

0.29 

158 

22 

4 

55 

0.1-0.15 

25-40 

2 

C 

600 

0.44 

165 

35 

4 

50 

0.1-0.15 

25-40 

1-2 

D 

150 

1.5 

82 

59 

3-4 

75-100 

0.2-0.25* 

50-65 

1-2 

E 

200 

4.4 

52 

no 

4-5 

160-200 

0.4-0.5* 

100-130 

1-2 

F 

200 

5.2 

56 

138 

4-5 

150-190 

0.4-0.5* 

100-130 

1-2 


C. Weak reduction of the n = 3 step 

In FigjS7|A, we present the colormap of the bin counts obtained on a 200 nm long 
junction (not presented in the main text), as a function of the dc and rf drives I and 
Rf. The measurement is taken at f = 3 GHz and 800 mK. For this frequency, we observe 
the complete suppression of the first Shapiro step. Again, the first oscillation is strongly 
disturbed as detailed in the main text. More surprisingly, a weak but distinct suppression 
of the third step n = 3 can be seen for low rf currents. It is emphasized by plotting 
the amplitude of each pair of steps (FigjST^). First, as described in the main text, the 
amplitude of the first step (red line) is fully suppressed up to the crossing point where 
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the n = 0 region vanishes, while the second step (green line) is fully visible. Then, the 
amplitude of the n = 3 step (cyan line) is smaller than the amplitude of the n = 4 one 
(purple line) while the opposite is usually expected from RSJ models. 

This could constitute evidence of a missing or suppressed third step, clearly observed 
only in one sample. However, other irregularities in the higher order steps are sometimes 
visible (see Fig.3 in the main article), so that this evidence has to be considered carefully. 
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FIG. S7. Weak reduction of the n = 3 step - (A) The differential resistance dV/dl is 
plotted as a colorscale, as a function of the dc and rf current drives / and I^f for a frequency 
of / = 3 GHz. A weak suppression of the third step (n = 3) is observed for low rf currents. (B) 
The step amplitudes are extracted by binning voltages (see main text) and presented for steps 
0 to 4 (with an offset of 0.03 for n = 1, 2 and 0.6 for n = 3,4). As commonly observed in our 
samples, the first step (red line) if fully suppressed before the crossing point. More interestingly, 
the amplitude of the third step (cyan line) is reduced as compared with our usual observations. 
In particular, it remains smaller than the amplitude of the fourth step (presented in magenta) 
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D. Hysteresis, bias instability and Shapiro steps 


In FigjSSl we present two datasets measured on the same junction as the one presented 
in section [TTC at two different temperatures. To emphasize the role of the control param¬ 
eters (dc and rf currents / and I^/), we now plot the differential resistance dV/dl as a 
colormap, as a function of these two parameters. On the left panel, measurements were 
obtained at / = 3 GHz at the base temperature of the dilution refrigerator (T = 12 mK). 
On the right panel, measurements were obtained in the same conditions (/ = 3 GHz) ex¬ 
cept for the temperature, here set high enough to suppress the hysteresis (T ~ SOOmK). 

Shapiro steps are identified as black regions where dV/dl ~ 0, while blue lines between 
black regions emphasize transitions between the different plateaus. A simultaneous read¬ 
ing of the voltage V gives access to the step index n (a few of them are indicated directly 


on Fig,S8). The dc current is swept in the direction indicated by the white arrow (from 
negative towards positive bias). At T = 12 mK, the very clear asymmetry at low rf excita¬ 
tions is a signature of the hysteresis observed in the I-V curve presented in the main text. 
The bistable dynamics that leads to hysteresis prevents the development of the phase- 
locked dynamics responsible for the Shapiro steps. Consequently, the latter are missing 
in the hysteretic region. In contrast, all steps are clearly visible at higher temperature 
(T ~ SOOmK), except for the n — 1 step which is fully suppressed ai f — 3GHz. 

It is possible to make a clear distinction between missing steps due to hysteresis and 
the missing n = 1 step attributed to a dvr-periodic supercurrent. First, as can be seen 


on Fig.SS, hysteretic switchings present vertical tangents (similar to ref l77|) while the 
doubled Shapiro step always exhibits a finite slope. Second, hysteresis is characterized by 
an asymmetry depending on the sweep directions, which is not the case of the missing step 
in our measurements. Third, the anomalous splitting of the n = 1/n = 2 steps beyond 
the pinch-off of the supercurrent (corresponding to the dark fringe in the bin counts, 
discussed in the main text) that accompanies a missing n = 1 step remains visible at all 
temperature regardless of the presence of hysteresis for low rf powers. For these reasons, 
it appears clear that one can safely neglect hysteresis as the origin of the missing n = 1 
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step. To avoid any problem, most measurements were performed at a temperature high 
enough to suppress the hysteresis (typically 450 to 800 mK). 


T ~ 12 mK 

_y 
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FIG. S8. Effect of hysteresis on Shapiro steps - The differential resistance dVjdl is plotted 
as a colorscale, as a function of the dc and rf current drives I and I^f for a frequency of / = 3 GHz, 
at the base temperature of the dilution refrigerator (T = 12 mK) and at T = 800 mK. The 
asymmetry in the figure signals a strong hysteretic behavior in the I-V curve (see main text, 
Fig.2). The wide white arrow (left panel) symbolizes the sweep direction of the dc bias current 
I. In the hysteretic region, low index n Shapiro steps are not visible. However, for high rf power, 
hysteresis vanishes and Shapiro steps reappear. The typical splitting of the n — 1/n — 2 steps 
is in particular still observable. In contrast, all steps (apart from the missing n — 1) are visible 
at T = 800 mK, and no asymmetry is observed. 
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E. Shapiro steps on a shunted device 


Furthermore, it is in fact possible to rule out definitively bias instabilities as a possible 
mechanism for the missing n = 1 step. By adding a shunt resistance in parallel with the 
junction, one can indeed suppress hysteresi^ and work in a configuration that approaches 
the voltage bias regime. To do so, we add a 10 Q resistor in series with the junction, and 
shunt these two elements with a 1 resistor. Thus, the current flowing through the 
junction is accessed by measuring the voltage across the 10resistor, together with the 
voltage across the junction. 



DC current I [pA] 
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FIG. S9. Shapiro steps on a shunted device - (A) I — V curves in the shunted configuration, 
at temperature T=12mK. A measurement in the absence of rf irradiation is shown as a black 
line. The low voltage regime with dV/dl < 0 corresponds to an unstable region when a pure 
current bias is applied. The colored lines present data measured under irradiation at / = 
4.2 GHz. The rf current is increased progressively from indigo to yellow lines. (B) and (C) 
Histogram of the voltage and differential conductance dl jdV as a function of DC voltage V and 
rf current Irf. As shown in the standard current bias configuration, a clear vanishing of the 
n = 1 Shapiro step is observed. 

Results obtained at the base temperature of the fridge (12 mK) and with an rf excita- 
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tion at 4.2 GHz are shown in Fig,^ First, the black solid line in FigjS9|A. represents the 
I — V curve in the absence of rf irradiation. At low voltages, one sees that dV/dl < 0. 
In this region, a pure current bias generates bias instability and hysteresis, which is here 
suppressed by the shunt resistance. When the rf irradiation is switched on (colored plain 
lines), Shapiro steps become visible in the I — V curves. As previously, the step n = 1 is 
clearly suppressed at low rf excitation (blue lines). Increasing the rf drive (from blue to 
yellow lines), one sees the n = 1 step is recovered at high drive amplitude as previously. 


Fig,S9 B and C present the voltage histograms and differential resistance dl/dV as a 
function of the normalized DC voltage V and rf current A/. Shapiro steps are visible as 
previously as maxima following horizontal lines. As for the measurements shown in the 
main text, the first step n = 1 is fully suppressed up to the oscillating regions. Though 
the contrast is not as good, the ’’dark fringe” at finite voltage described in the main text 
is also visible. 


F. Magnetic field dependence 

Further investigations of the anomalous Shapiro response have been carried out in 
the presence of perpendicular-to-plane magnetic fields. First, when the I-V curve is 
measured without rf irradiation, a Fraunhofer-like diffraction pattern of the critical current 
is observed (plotted as a red line in Fig jS10| upper panel). In the junction presented here, 
the pattern is slightly distorted (probably due to flux-trapping in the magnet). The 
periodicity in the magnetic field has been evaluated from undistorted patterns in various 
other samples. It corresponds to a conventional periodicity as previously reported on 
similar sampled, and as expected for ballistic systems with such aspect ratioP^^^. 

The observation of Shapiro steps in the presence of the magnetic field reveals that the 
anomalous doubling of the first step is only observed when the critical current is high 
enough, namely only at the center of the central lobe (B < 0.3 mT) and to a lesser extent 
at the centers of the second and third lobe. The ratio Q 12 of the first (n = 1) and second 
(n = 2) step widths (as introduced in the main text) is plotted as a function of the 
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FIG. SIO. Magnetic field dependence - (upper panel) Critical current Ic as a function of the 
perpendicular magnetic field B. The typical diffraction pattern of the supercurrent is observed, 
(lower panel) The ratio Q 12 is presented as a blue line, as a function of the perpendicular 
magnetic field B. It follows the same oscillations as the critical current Ic- The grey area 
between the dashed lines show the region in which the ratio Q 12 falls, as evaluated by a RSJ 
model. 


magnetic field in Fig,SIO (as a blue line, lower panel). The plot reproduces the variations 
observed in the critical current: the reduction of the first step is distinctly observed for 
high critical currents, yielding low values for Q 12 , while Q 12 increases for low critical 
currents near the minima of the Fraunhofer pattern. This behavior can be understood 
by assuming a constant value of the dyr-supercurrent fraction Ia-k/Ic'- fhe frequency Bn 
rapidly decreases when B decreases, so that a correct observation of the anomalous step 
is not possible according to the criterion / < f^n- 

One could expect that a small magnetic field induces a splitting of the Andreev bound 
states. This would consequently weaken the effect of Landau-Zener transitions (see section 


IV C) and allow us to exclude or confirm Landau-Zener transitions as responsible for the 


anomalous Shapiro response. Given the large p-factor in HgTe (around 5 - ~ 20 for bulk 
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HgTe), this effect could show up at relatively weak helds. To generate a splitting of 0.01 Aj 
that significantly alters Landau-Zener transitions, a typical magnetic field Bq is required, 
with qubBo ~ O.OlAj, yielding Bq ~ 3.5 mT. The presence of a week magnetic field 
does not seem to suppress the dvr-periodic contribution, ft thus does not reveal any clear 
signature of energy splitting of the Andreev spectrum. This supports an interpretation 
based on the presence of a topological Andreev bound state. Indeed the latter is predicted 
to have no spin-degeneracy and therefore the level crossing at phase differences vr, Svr, etc. 
should persist even in the presence of a magnetic field on the TI. 


III. ADDITIONAL RESULTS ON GRAPHENE-BASED JUNCTIONS 

A. Typical properties of the graphene-based junctions 

In this section, we briefly present measurements of the Shapiro response of graphene- 
based devices. First, graphene flakes on hexagonal Boron-Nitride are fabricated by a van 
der Waals stacking method. Superconducting contacts are patterned by standard electron 
beam lithography and Niobium is deposited by magnetron sputtering. The geometry of 
the presented devices is similar to the one used for the HgTe devices, with the graphene 
weak link width ranging between 1.5 pm and 2.5 pm, and the length between 200 nm and 
300 nm . 

The graphene flake exhibit densities between —2.5 x 10^^cm“^ and 2.5 x lO^^cm”^, 
which can be tuned by the means of a back-gate. In the left panel of Fig jSll[ we present 
the extracted normal state resistance Rn of the device as a function of the back-gate 
voltage Vg. The mobility approaches 3000 - 8000 cm^ s“h Close to the Dirac point, 
this corresponds to a mean free path of around 50 nm. Consequently, these devices are 
not in the ballistic limit but are relatively close to it. 

In the right panel of Fig jSll| we present the differential resistance dV/dl of one junction 
as a colormap, as a function of the back-gate voltage Vg and dc drive current I. The 
supercurrent then appears as a dark region (dV/dl ~ 0), while the normal state exhibits 
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FIG. Sll. DC characterization of Josephson junctions with graphene weak link - 

(left) Normal state resistance of a graphene-based junction as a function of back-voltage Vg. 
The Dirac point is clearly indicated by a maximum around Vg ~ —IV. (right) Colormap of 
the differential resistance dV/dl as a function of the back-gate voltage Vg and dc drive current 
I. The supercurrent is visible as a black region (dV/dl ~ 0) over the whole back-gate range, 
ranging from a few tens of nA around the Dirac point up to 1.8 pA for large gate voltages. 


finite values of dV/dl. The devices show a supercurrent over the whole back-gate range, 
yielding supercurrents of a few tens of uA around the Dirac point and up to 1.8 pA for 
large gate voltages. 


B. Shapiro response of the graphene-based junctions 


In this section we present two typical sets of data taken on the graphene-based junctions 


(see Fig,S 12 ). The differential conductance dl/dV is plotted as a colorscale, as a function 
of the dc voltage V and rf current drive 7^/ for two frequencies (/ = 7GHz and 5.5 GHz). 
Despite fluctuations in the measurements (originating in fluctuations of the offset of one 
of our amplifiers), clear Shapiro steps are visible at the expected voltages 14 = nhf /2e. 
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As frequency is lowered, the steps become less discernible but no sign of a vanishing n = 1 
step is observed in these sets. 

Measurements have been carried out for different values of the back-gate voltage (Dirac 
point, high n or p doping) and down to 4 GHz (for which the steps are hardly visible) but 
no sign of a missing n = 1 step has ever been spotted. Since no dyr-periodic modes are 
expected in the graphene junction, the crossover freqnency is expected to be / 4 ^ = 0 in this 
system. These measurements establish an experimental upper boundary / 47 r <C 4 GHz. 
Given the characteristic Josephson frequency of the device fj ~ 72 GHz, this sets an 
upper boundary of <C 75 nA. 
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FIG. S12. Shapiro response in graphene-based junctions - The differential conductance 
dljdV is plotted as a colorscale, as a function of the dc voltage V and rf current drives Irf 
for two different frequencies / = 7 GHz and 5.5 GHz. The data has been taken at the base 
temperature of the dilution refrigerator (T = 25 mK), close to the Dirac point. All Shapiro 
steps are visible, up to |n| > 12. 
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IV. ELEMENTS OF MODELING 


A. Simulations using RSJ equations 


Simulations based on a standard RSJ modeP^ have been carried out to compare our 
results with a simple and well-understood model. In these simulations the junction is 
modeled together with a resistive shunt to capture the impedance of the environment 
(that plays an essential role in the dynamics of the junction). In this framework, the total 
current through the system / can be written as the sum I = Ir + Is where Ir = ^ 
is the current through the resistor and Ig the supercurrent through the junction!^. 
Combining the hrst Josephson equation dcjy/dt = 2eV/h, and the current bias I = Idc + 
Rf sm2nft obtains a first order ordinary differential equation: 

—^0 Is{(p) = Idc + Irf sin27r/t (2) 

ZCrCn 


The rf excitation is here represented as a current Rf instead of a voltage V)./ in agreement 
with most of the literature on Shapiro steps. It assumes that the characteristic field 
impedance of the radiation field is high compared with the junction impedanc^. In our 
case (low capacitance), the typical impedance of the junction is given by its resistance 
(typically between 30 to 150 fl), smaller than the free space impedance (in which RF 
propagates before the junction) given by c/tq = I/ccq — 376 Ohm. This approximation is 
probably a bit crude, but it is further justified by the overall agreement of the effect of 
frequency on the Shapiro response (Fig. 3 of the main text and Fig. S5 and Fig. S17[ ), 
which is a characteristic feature of the current bias model developed by Russei®'. We 
simulate the results of this equation using a simple Runge-Kutta algorithm (RK4) to 
obtain the I-V curve, realize a binning of the voltage V and finally compute the ratios 
Qi 2 ;Q 34 for comparison with experimental results. 

The results shown here are obtained with a 27r-periodic sinusoidal current-phase rela¬ 
tion Is{(f>) = Rsinf. The results obtained by Russei^^are reproduced. The global visual 
agreement with our measurements is good, except the missing first step not accounted for 
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FIG. S13. Simulated Shapiro response of a conventional Josephson junction - In the 

framework of the RSJ model (Eqj^, we obtain colormaps of the bin counts as a function of the 
voltage bins (in normalized units) and rf current drive I^f. For a high frequency (/ = 0.4/j) 
the results obtained resembles that of Fig.3C in the main text. All steps are clearly visible and 
appear one by one as A/ increases. As frequency / is decreased, the pinch-off of the supercurrent 
is moved to higher powers. One also observes a decreasing period in the oscillating pattern. 


in the standard RSJ model. In particular, the increasing frequency of oscillations (with 
respect to the rf current needed for the pinch-off of the supercurrent) is clearly observed. 

In the lower panels, the step amplitudes for n = 0,1, 2, 3,4 are plotted as a function of 
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the rf current amplitude Irj- Again, one observes the effect of the excitation frequency /. 
As / decreases, the width of the first lobe gets much larger than the widths of the other 
lobes. Moreover, these graphs also show that the step amplitudes (hence their visibility) 
decrease (with respect to the critical current). This phenomenon, clearly observable in our 
measurements (lower panels of Fig.3, main text), limits our measurements to / > 2 GHz. 


B. Extended RSJ model with a 47r-periodic contribution: 


Simulations have also been carried out when adding a dyr-periodic contribution follow¬ 
ing Dominguez et al^. We present in this section our results : the appearance of the 
doubled step is qualitatively well described as in the previous reference, but a quantitative 
agreement has not been obtained yet. The results are presented below in the following 
manner. First, we show that the addition of a dyr-periodic contribution / 4 ^ sin 0/2 to a 
sinusoidal 27r-periodic current-phase relation (CPR) l 2 -„ sin 0 is responsible for the dis¬ 
appearance of all odd steps at low frequency / < in a comparable way to what we 
experimentally observe on the first step n = 1. Then the marginal effect of the 27r-periodic 
CPR is illustrated by comparing I-V curves in the presence of rf irradiation for a few dif¬ 
ferent CPRs. No generic 27r-periodic CPR is fonnd to show missing odd steps, while the 
addition of a small dvr-periodic term enforces the disappearance of odd steps regardless 
of the 27r-periodic term. 

a. Effect of frequency on a 27 t + Air-periodic supercurrent: Using Eq.([^, we simulate 


the effect of a small dvr-periodic contribution sin 0/2. In Fig.Sld, we show the effect of 
frequency on a CPR of the form sin (l)/2 + E-k sin 0 by plotting for different frequencies 
the histograms of the voltage as fnnction of the rf current Ef (upper panels) and the 
amplitndes of steps n = 0 to n = 4 (lower panels). The parameters are Et^ = 1, Ett = 0.15, 
so that 04 ^ = 0.15/j. For a high freqnency / > /47r (left, / = 0.5/j), all steps are visible, 
and the 2D plot is similar to the one obtained without 4vr-periodic contribution. As 
frequency is lowered (center, / = /47r = 0.15/j, and right, / = 0.05/j), the amplitude of 
odd steps progressively decreases, and these odd steps are completely suppressed at low 
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power. The simulated behavior of the n = 1 step is very similar to the one experimentally 
observed, and thus reinforces our interpretation. Though the crossover is not very sharp, 


it happens for / = / 4 ^ (= 0.15/j in the case of Fig,S14). One also observes that the 
oscillatory pattern (at high power) is also progressively modified from a 27r- to a dyr- 
dominated pattern. In particular, odd steps show a very pronounced first minimum. The 
dark fringe we experimentally observe in the oscillatory pattern is understood as the result 
of the progressive towards a pattern with a halved number of oscillations, thus yielding 
progressively suppressed lobes. 

b. Effect of the CPR Though the presence of a small dvr-periodic contribution 
/ 47 rSin 0/2 is found necessary to observe vanishing odd steps, the exact description of the 
27r-periodic supercurrent does not infiuence much the Shapiro response. To illustrate this 
finding, we focus on three different CPR. The first one is a standard R-k sin f contribution 
(obtained for a tunnel junction for example), as in the previous paragraph. For the other 
two, we assume that the current is carried as a single mode of transmission r = 0.8 and 
0.99 respectively. The skewness of the CPR increases with transmission. The current 
carried by such a mode is then given by rsin^/y/l — r sin^ <p/2. We normalize all CPRs 


such that /c = / 27 r = 1. Fig jS15| shows in the upper panels the CPR and in the lower pan¬ 
els an I-V curve in the presence of rf irradiation, for a low frequency / = 0.05/j = / 47 r/ 3 . 
On the left side, the CPR contains only this 27r-periodic supercurrent, and the obtained 
I-V curves exhibits all integer steps and do not depend on the exact CPR. When a 
dvr-periodic supercurrent is added, all odd steps are suppressed or reduced. The CPR 
matters marginally, and only seem to control the exact detail of the transition from the 
2 vr- to the dvr-dominated regime : the transition is slightly faster for a skewed CPR. 


C. Evaluation of Landau-Zener probabilities 


If some Andreev bound states have very high transmission (or small energy gap 5, 


see Fig.Slfi), non-adiabatic Landau-Zener processes may be responsible for an additional 


dvr-periodic contribution to the supercurrent. Namely, some gapped 2 vr-periodic Andreev 
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FIG. S14. Simulated Shapiro response with a 47r-periodic term - (upper panels) Col- 
ormaps of the bin counts as a function of the voltage (in normalized units) and rf current drive 
Irf^ in the presence of a small dvr-periodic contribution (such that / 47 T = 0.15/j). (lower panels) 
Amplitudes of step n = 0 to n = 4, as a function of the rf current drive Aj. At high frequency 
(/ = 0.5/j > / 47 r), all steps are visible, and the result is very similar to the one obtained without 
dvr-periodic modes. As the frequency is decreased (/ = = 0.15/j and / = 0.05/j), the am¬ 

plitude of the odd steps decreases, in particular the n — 1 step which is completely suppressed 
at low power. Simultaneously, the oscillatory pattern (high power) develops anomalies, with 
very pronounced first minima on the odd steps. 


levels would behave as effective 4vr-periodic modes, in the absence of true 4vr-periodic 
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FIG. S15. Effect of the CPR - (upper panels) Investigated current-phase relations: in blue, 
the sinusoidal 2-periodic supercurrent I 2 tt sin 0, in red and green, single-mode-like supercurrent 
with transmission r = 0.8 and 0.99 respectively. On the right side, an additional dvr-periodic 
contribution is added, with amplitude 1/^^^ = 0.15/27r. (lower panels) I-V curves in the presence 
of rf excitation. While all steps are visible in the absence of 4-periodic modes, odd steps are 
missing or strongly suppressed when they are present. The exact CPR plays a small role in the 
crossover from the dvr-dominated regime to the conventional one when frequency is increased. 

modes. To our knowledge, studies of the effect of Landau-Zener transitions on the dy¬ 
namics of Josephson junctions are scarcd^^^^^^ in general and systems with multiple An¬ 
dreev levels have not been explored yet. We thus use a single mode approximatioipS^ and 
crudely assume that one level (with lowest gap 5) has a predominant role in Landau-Zener 
processes. 

Following Dominguez et a/.^, we include in our simulations stochastic Landau-Zener 
transitions at the anticrossings, with a probability P, and partially reproduce their results. 


dl 













The exact motion of the phase difference is hard to picture, but the general trend can be 
understood on a heuristic basis. The phase will undergo shifts of 4?! (per period of the 
drive) in the case of a Landau-Zener transition and only 27r in the absence of transitions. 
On average this results in a non-universal Shapiro step that is neither hf /2e nor hf/e. 
These two values are recovered for P = 0 and P = 1 respectively. The results of our 
simulations are presented in Fig jS17| A-, in which a close-up of the n = 2 Shapiro step is 
presented, for different values of P. As Dominguez et ai, we observe a splitting of the 
step for P close to 1, with voltage plateaus that deviate from the quantized value hf/e 
until the split steps are eventually less discernible for P < 0.7. Our experimental results 
do not show any splitting or deviation to the quantized value, with an accuracy of a few 
percents. This indicates that the probability must be P > 0.97 (P = 1 being equivalent 
to having a fully gapless mode). 



FIG. S16. Effect of Landau-Zener transitions - A typical spectrum of gapped 27r-periodic 
Andreev levels is presented as a function of the superconducting phase difference f. In the 
presence of Landau-Zener transitions at the anticrossings (represented as green arrows for f = 
tTjStt, ...), some energy levels with small gap 5 (high transmission) could contribute as a 47r- 
periodic component in the supercurrent. 
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From the previous model, one can estimate the importance of the Landau-Zener tran¬ 
sitions. First, we solve Eqj^ and obtain the phase and its derivative <p{t) oc V{t) 
as a function of time t. The parameters are chosen such that the junction lies on the 
first Shapiro step (V = {V(t)) = |^). In that case, a high probability of Landau-Zener 
transitions would lead the junction to exhibit a doubled step. 

A typical plot of (j){t) and V{t) is shown in Fig jSlTp , with (j){t) as a red line and 
V(t) oc 4>{t) as a blue line for the following parameters: / = 0.5 Ic, hf = 0.8 /c, / = 0.2 fj. 
One first observes that the phase 0 follows an anharmonic motion synchronized with the 
excitation drive at frequency /: during one period of duration 1//, the phase (f) increases 
by 27r, yielding an averaged voltage V = hf/2e as expected for the first Shapiro step. 
Equivalently, one can calculate the average of V{t) and obtain V = {V{t)) = |^. Then, 
we access the time t for which 0 reaches the anticrossing (for (pit) = Stt for example) 
and read the derivative of the phase ^Is^r at this point or equivalently the voltage V\z-w 
Finally, the Landau-Zener probability P is obtained from the following equatiorpSl 

P = exp ( — 27r-r 

This assumes a generic pair of Andreev bound state with energy e±{(j)) = + A? cos^ 0/2. 

For the graphs presented here, we obtain V\ 2 ,-k — 6-4 |y. This yields a probability P > 0.98 
if 5 < 6peV. The most stringent constraint is obtained for the lowest frequencies acces¬ 
sible in the experiment (2GHz). With the criterion P > 0.97, we obtain S < 9peV for 
Aj = 0.12 meV (L = 400 and 600 nm), and S < 18peV for Aj = 0.35 meV (L = 150 nm). 
This corresponds to a transmission >/l — ((j/Aj)^ > 0.994 in both cases. 
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FIG. S17. Landau-Zener transitions in the RSJ model - A) Evolution of the Shapiro step 
n = 2: I-V curves obtained for different values of the Landau-Zener probability P are plotted, 
with a focus on the n — 2 Shapiro step. When P = 1, the second step (n = 2) is fully developed 
and reaches the expected quantized value hf /e. For P < 1, it progressively deteriorates: the 
voltage departs from the quantized value /i//2e, and the plateau becomes less visible. B) Time 
evolution of phase difference and voltage V{t)\ The phase difference <^{t) (red line, in units 
of tt) and voltage V{t) (blue line, in units hf/2e) across the junction is calculated via RSJ 
equations and plotted as a function of time t in units of the rf period 1//. The estimation of 
the voltage EIstt at = Stt (anticrossing) enables a proper evaluation of the Landau-Zener 
transition probability. 
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